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^! Abstract 

The Overlap-Dirac operator provides a lattice regularization of massless vector 
' gauge theories with an exact chiral symmetry. Practical implementations of this 

. operator and recent results in quenched QCD using this Overlap-Dirac operator are 

reviewed. 
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Since tlie beginning of lattice gauge theory the regularization of chiral fermions 
^ ' has been afflicted with severe problems. When regulating fermions on a lattice, 

p • typically, unwanted doublers with opposite chirality appear. These doublers 

O . can be lifted (given mass of the order of the cut-off) at the cost of explicitly 

^ ' breaking chiral symmetry, as in the Wilson fermion formulation. Alternatively, 

a remnant of chiral symmetry can be retained, with a smaller number of dou- 
^ ' biers interpreted as flavors, as in the staggered fermion formalism. However, at 

■ flnite lattice spacing the flavor symmetry is broken. Both these approaches fail 

from the outset for regulating Weyl fermions. The central problem in the non- 
perturbative regularization of gauge theories with Weyl fermions is to write 
down a formula for the fermionic determinant when the fermion is in some 
complex representation of the gauge group, since depending on the topology of 
the gauge field the chiral Dirac operator can be square or rectangular, where 
the difference between rows and columns is the index. 



Significant progress in the formulation of chiral gauge theories has been made 
by the overlap formalism [1]. The overlap formalism was inspired by two pa- 
pers [2]. The central idea is that an infinite number of Dirac fermions (la- 
beled by s) with a mass term of the form iPs{PlMss' + PrMI^,)iPs' and chiral 
projectors Pl,r can be used to regulate a single Weyl fermion if the infi- 
nite dimensional mass matrix M has a single zero mode but has no zero 
modes. Kaplan's paper [2] uses this idea to put chiral fermions on the lattice 
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where they are referred to as domain wall fermions since Kaplan used a mass 
matrix that has a domain wall like structure. In the overlap formalism the infi- 
nite number of Dirac fermions is described by two non-interacting many body 
Hamiltonians, one for each side of the domain wall, and the chiral determinant 
is written as the overlap between their groundstates 



|0— ) is the many body ground state of 7i~ = a^-^^a and |0+) the many body 
ground state of Ti"*" = Hyj{U)a, with ■y^H^iU) = D^{U) the usual Wilson- 
Dirac operator on the lattice with a fermion mass in the supercritical region 
{rric < m < 2). a (a^) are canonical fermion annihilation (creation) opera- 
tors. On a finite lattice, the single particle Hamiltonians are finite matrices 
of size 2K x 2K with K = V x N x S where V is the volume of the lat- 
tice, is the size of the particular representation of the gauge group and 5* 
is the number of components of a Weyl spinor. Then |0— ) is made up of K 
particles. If |0+) is also made up of K particles, then the overlap is not zero 
in the generic case. If the background gauge field is such that there are only 
K — Q negative energy states for Hw{U) then the overlap is zero. Any small 
perturbation of the gauge field will not alter this situation. Furthermore, the 
overlap (0 — |aj^ ■ • ■ 0-1^10-1-) will not be zero in the generic case if the fermion 
is in the fundamental representation of the gauge group showing that there is 
a violation of fermion number by Q units. So, clearly, the overlap definition 
of the chiral determinant (1) has the desired properties. 

A generic problem with simulations of chiral gauge theories is that the chiral 
fermion determinant is complex. A "brute force approach" is feasible in the 
simulation of two dimensional models, and in this way the overlap formal- 
ism has successfully reproduced non-trivial results in two dimensional chiral 
models on the lattice [3]. The brute force approach, however, is clearly not 
feasible in four dimensions, where efficient numerical techniques are essential. 
This prevented simulations of chiral gauge theories and tests of the overlap 
formalism in four dimensions so far. 

Clearly, any formulation of lattice chiral gauge theories is also a formulation 
of massless vector gauge theories with an exact chiral symmetry and a pos- 
itive fermion determinant (the product of the chiral determinant for the left 
handed fermions and its complex conjugate for the right handed ones). Lattice 
QCD using the overlap formalism reproduces the well-known mass inequali- 
ties between mesons and baryons, and the U{Nf)v x U{Nf)A symmetry in an 
Nf flavor theory is broken down to U{l)v x SU{Nf)v x SU{Nf)A by gauge 
flelds that carry topological charge [4]Q. If the SU{Nf)v x SU{Nf)A symme- 
try is spontaneously broken, then massless Goldstone bosons should naturally 

^ See section 9 of [1] for details. 



det C{U) ^ (0 



0+). 




2 



emerge in the overlap formalism. Since the symmetry breaking pattern is ex- 
actly as in the continuum, all the soft pion theorems should hold on the lattice 
as well. 



For a vector gauge theory the computation of the fermionic determinant can be 
simplified significantly compared to the original version (1) that involves the 
computation of the overlap of two many body ground states. One way to derive 
the simplified expression is to start with the variant of domain wall fermions 
of ref . [5] , apphcable for a vector gauge theory. We choose this approach here 
to emphasize the close connection between domain wall and overlap fermions. 
Integrating out all the fermion and Pauli-Villars fields, Ncuberger derived 
the following expression for the determinant describing a single light Dirac 
fermion [6] 



det Dowil^] Lg) — det 



1 + /i + (1 — /i)75 tanh 



In T,„ 



(2) 



Here is the transfer matrix in the extra direction, whose extent, Ls, has 
been kept finite, and < /i < 1 describes fermions with positive mass all the 
way from zero to infinity. In [5] the fermion mass /i is denoted by m/. In the 
limit Lg ^ oo (2) becomes 

detr>w(//) =det|^[l + /x + (l-/i)75e(-lnr^)]|. (3) 



It is only in this limit that massless domain wall fermions have an exact chiral 
symmetry. Finally, taking the lattice spacing, a^, in the extra direction to zero 
one obtains the Overlap-Dirac operator of Neuberger [7] 

D{,i) = ^[l + fx+{l-fx)j,e{H^)]. (4) 



The external fermion propagator is given by 

D-\^) = {1 - ^)-' \d-\i,) - l] . (5) 



The subtraction at = is evident from the original overlap formalism [1] 
and the massless propagator anti-commutes with 75 [4,7]. With our choice of 
subtraction and overall normalization the propagator satisfies the relation 

pi{b^ [^,D-\pi)f\b) = {b^D-\pi)\b) V b satisfying ^,\b) = ±\b) (6) 



for all values of /i in an arbitrary gauge field background [8]. The fermion 



3 



propagator on the lattice is related to the continuum propagator for small 
momenta and small /i by 

D;\m,) = Z^'D-\^) with m, = Z^V (7) 

where and are the mass and wavefunction renormalizations, respec- 
tively. Requiring that (6) hold in the continuum results in Z^Z^n = 1. We find 
that a tree level tadpole improved estimate gives 

Z^^Z-'^-[m-A{l-Uo)] , (8) 

Uq 

where Uq is one's favorite choice for the tadpole link value. Most consistently, 
for the above relation, it is obtained from mc, the critical mass of usual Wilson 
fermion spectroscopy. 

In the rest of this paper we discuss practical implementations of the Overlap- 
Dirac operator and present some recent results. Owing to the recent flurry [9] 
of theoretical activity arising from the "unearthing" of the Ginsparg- Wilson 
relation [10], a few remarks are in order. The massless Overlap-Dirac operator 
in (4) satisfies the Ginsparg- Wilson relation 

^(0)75 + 75^(0) = 2Z}(0)75D(0) (9) 

implying that the massless propagator (^"^(O) — 1) anticommutes with 75. If 
we write 

(10) 

then the Ginsparg- Wilson relation reduces to = 1. Since we would want 
75D(0) to be Hcrmitian. /f„ shoiild be a Hcrmitian operator. With this reduc- 
tion of the Ginsparg- Wilson relation, it is easy to show that [11] 

detL'(O) = 1(0- |0+)f (11) 

i.e., the overlap formula for a vector theory. Here |0-|-) is the many body 
ground state of a^Had. This establishes a one-to-one correspondence between 
the overlap formula and the determinant of a fermionic operator satisfying the 
Ginsparg- Wilson relation for massless vector gauge theories. 



L>(0) = - l + ^^Ha 
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2 Practical Implementations of e{Hyj) 



In order to compute the action of e{Hyj) 



on a vector one can proceed in 



several different ways. Since we are interested in working in four dimensions, it 
is not practical to store the whole matrix H^. Therefore standard techniques 
to deal with the square root of a matrix [12] will not be discussed. 

One could attempt to solve the equation H^jCj) — Hyjh to obtain = e{Hyj)b 
using iterative techniques. Such techniques have been developed to solve linear 
systems with fractional powers of a positive definite operator using Gegenbauer 
polynomials [13] and applied to the Overlap-Dirac operator in [14]. 

Another approach is to efficiently approximate e{Hy,) as a sum of poles: 

N 



k=l ^ "fc 



(12) 



The action of e(if^) on a vector then involves a single conjugate gradient with 
multiple shifts [15]. 

One approximation, called the polar decomposition [16], has been adapted in 
this context and first used in the study of the three dimensional Overlap-Dirac 
operator by Neuberger [17]. Here the coefficient cq = and 

1 71 

Cfe = , ^ , , -; 4 = tan2 (2k - 1) . (13) 



In this approximation, 

(1 + zW + (l-z 



(l + ^)2iV_(i_^)27V 

e(z) ^ g^{z) = ) J J . (14) 



Clearly gN{z) = g^i^-lz) and 5'Ar(l) = 1. The error t{z) — gN{z) is strictly 
positive and monotonically decreases from z — Q to z — 1. 

For another approximation, called the optimal rational approximation [18], 
the coefficients are obtained numerically by an optimal fit using the Remez 
algorithm. The coefficients in a slightly different notation have been tabulated 
for A/" = 6, 8, 10 in Ref. [14]. We have found it necessary to use A^ = 14 for 
our recent applications. The coefficients in this case are shown in Table 1. In 
the optimal rational approximation, the approximation to t{z) has oscillations 
and is not bounded by unity. A plot of the approximation gN^z) obtained as 
a fit over the region [0.001, 1] is shown in Fig. 1 for A" = 6 to 14. While we fit 
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(2.413330975e+00,1.361747338e+01) 

(2.925707925e-01,1.213113539e+00) 

(1.166359792e-01,2.752627333e-01) 

(6.216038074e-02,6.543005714e-02) 

(3.423610040e-02,1.164228894e-02) 

(1.545550091e-02,9.937321442e-04) 

(3.419245947e-03,1.585925699e-05) 



(6.257184884e-01,3.135687028e+00) 
(1.737405612e-01,5.596349298e-01) 
(8.372555094e-02,1.364115846c-01) 
(4.652496186e-02,2.923946484e-02) 
(2.404754621e-02,3.887745892e-03) 
(8.436481876e-03,1.684882417e-04) 
(1.138166539e-03,5.914114023e-07) 



Table 1 

The coefficient pairs (c^, d^) k = 1,. . . ,N for the N = lA optimal rational approx- 
imation for which cq = 0.0850910. The coefficients are obtained as a result of an 
optimal fit over the interval [0.001, 1]. 
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Fig. 1. Plots of the optimal rational function approximation to e{z) (Eq. 12) over 
the interval [0.001, 1] for A'^ = 6, 8, 10, 12, and 14. The small z region of = 6 and 
8 is not shown. The fits are still good for some z > 1. 

over this region, the approximation is still good for z somewhat larger than 1. 
The approximation is bounded by unity only if 0.025 < z < 1.918 for N — 14. 
In this range the maximum deviation from unity is equal to 3.1 x 10~^. This 
range will increase if one increases the order of the approximation. For the 
current applications we found this range to be sufficient. 

The approximation to e{H^) by poles involves a multi-shift inner conjugate 

gradient and therefore it seems necessary to store N vectors where is the 
order of the approximation. One can avoid storing the extra vectors if one is 
willing to perform two passes of the inner conjugate gradient [19]. A Lanczos 
based algorithm that also avoids this extra storage by requiring two passes 
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has been proposed, but it involves an explicit diagonalization of a tridiago- 
nal matrix [20]. In this method is approximated by a small dimensional 
tridiagonal matrix (anywhere between 100 and 1000) and e{H^) is computed 
by first diagonalizing the tridiagonal matrix and then performing the trivial 
operation of e on the eigenvalues. The accuracy is increased by increasing the 
order of the tridiagonal matrix. 

Since in practice it is the action of D{ij,) on a vector we need, we can check for 
the convergence of the complete operator at each inner iteration of e{Hyj). This 
saves some small amount of work at /i — and more and more as /i increases, 
while at /i = 1 (corresponding to infinitely heavy fermions, c.f. Eq.(4)) no 
work at all is required. 

Each action of e{H^) involves several applications of Hy^ on a vector with the 
number depending on the condition number of Hy^{m) in the supercritical mass 
region. In Figure 2 we show the density of (near) zero eigenvalues for m = 1.7. 
We see that while p(0; 1.7) decreases rapidly as P increases, it does not appear 
to go zero at a finite lattice spacing. The second part of Figure 2 emphasizes 
that p(0; 1.7) decreases exponentially in some power of 1/a (the power here 
is not well determined, but reasonably fits 1/2). This result imphes that on 
a specific gauge background, H^j could have an arbitrarily large condition 
number due to a few small eigenvalues. This can make the computation of 
(f) = e{Hyj)h expensive for all methods considered. In addition some care is 
needed when using the approximation of e(if^) by a sum over poles. Clearly 
e{Hy,) can be replaced by e{sHy,) where s > is an arbitrary scale factor. 
Wc should choose the scale so that the maximum eigenvalue of sH^ is not 
above the range where the approximation is deemed good. Having so chosen a 
value for s, wc need to deal with the low lying eigenvalues that fall outside the 
range of the approximation to e(i/^). We do this by computing a few low lying 
eigenvalues and eigenvectors of if^, for which we then know the contribution to 
e{Hy^) exactly, and projecting them out before applying the approximation to 
the orthogonal subspace for which the approximation is good. The number of 
eigenvalues that have to be projected out will depend on the lattice coupling, 
the lattice size and the lower end of the range of the approximation. It will 
roughly increase with the volume at a fixed coupling making it difficult to go to 
large lattice volumes at strong coupling. However, the number of eigenvalues 
that have to be projected out will decrease as one goes to weaker coupling 
even at a fixed physical volume. This is because the density of eigenvalues of 
near zero goes to zero as one goes to the continuum limit [21]. The Ritz 
functional method [22] can be used to efficiently compute the necessary low 
lying eigenvalues and eigenvectors of H^j. 

For the case of domain wall fermions at finite extent Lg in the fifth direc- 
tion, the degree to which tanh(— Ls(lnTu,(m))/2) approximates e(— InT^(m)) 
is determined by Lg and the eigenvalues of Tyj{m) near 1. One can show ana- 
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Fig. 2. The approach of p{0; 1.7) to the continuum hmit versus j3 and versus l/vc?a. 

lytically [1] that in a fixed gauge background a unit eigenvalue of T^{m) and 
a zero eigenvalue of H^{m) occur at the same mass m. Also, the change of the 
corresponding eigenvalue in m is the same for both T^{m) and H^{m). This 
imphes that the density of zero eigenvalues p(0; m) is the same for both Hyj{m) 
and — ln(T^(m)). The degree to which these zero eigenvalues affect physical 
results is determined by the physical observable, L^, and the fermion mass 
Studies of the Lg dependence for various quantities at non-zero fermion masses 
can be found in Ref. [23]. In particular, larger Lg at fixed mass n is needed 
for stronger coupling due to the increasing p(0; m). 



3 Spontaneous chiral symmetry breaking 



Spontaneous chiral symmetry breaking is an important feature of QCD. How- 
ever, it is not fully realized on a finite lattice and at finite quark mass. Thus 
one needs to carefully study the approach to the infinite volume and chi- 
ral limit. Conventional lattice fermion formulations explicitly break the chi- 
ral symmetry (at least partially) at finite lattice spacing, obscuring the ap- 
proach to the infinite volume and chiral limit in practical simulations. Overlap 
fermions preserve the chiral symmetry at finite lattice spacing. This should 
facilitate a study of spontaneous chiral symmetry breaking. As a practical test 
of the Overlap-Dirac operator, we consider quenched QCD. The chiral limit 
of quenched QCD is tricky, though, because topologically non-trivial gauge 
fields are not suppressed in this limit. Gauge field topology results in exact 
zero modes of -D(O) as long as one is in the supercritical region of if^(m). This 
is demonstrated in in Fig. 3 where we show the spectral flow of eigenvalues of 
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Fig. 3. (a) The low lying spectrum of Hg (octagons) in a pure SU(2) gauge field 
background at /3 = 2.5 together with the spectral flow of Hyj (lines). The zero mode 
of Ho found for m > 0.8 is singly degenerate and is associated with an instanton in 
the background gauge field, (b) Shape of the zero mode of Ho projected to one axis 
for the various m (octagons) and of H^ at the crossing point (line) . 

both Hw{rn) and Ho = j^Di^O) as a function of m for an SU(2) gauge back- 
ground at /3 = 2.5 on an 8^ lattice. We see a single level crossing zero near 
m — 0.9 in the spectral flow of H^,. At this mass, we see the sudden appearance 
of a zero eigenvalue (with chirality 1) among the smoothly changing non-zero 
eigenvalues (in opposite sign pairs with chirality equal to their eigenvalue). 

Zero eigenvalues of Ho due to global topology have a definite chirahty. The 
spectrum of Ho is in [-1,1] and the non-zero eigenvalues of Ho that have a 
magnitude less than one come in pairs, ±A. The associated eigenvectors are not 
eigenvectors of 75, but rather 75 has expectation value ±A in the eigenvectors, 
ip'^'j^ip = ±A. Since Ho is an even dimensional matrix, the unpaired zero 
eigenvalues have to be matched by unpaired eigenvalues equal to ±1. This is 
what is expected to happen in a topologically non-trivial background. It is 
straightforward to obtain the spectrum of D(n) from the spectrum of Ho{0). 
Due to the continuum like spectrum one can study the approach to the chiral 
limit using the Overlap-Dirac operator by separating modes due to global 
topology from the remaining non-zero eigenvalues [8]. 

The main quantity that needs to be computed numerically is the fermion prop- 
agator D~^{fi) in Eqn. (5). Certain properties of the Overlap-Dirac operator 
enable us to compute the propagator for several fermion masses at one time 
using the multiple Krylov space solver [15] and also go directly to the massless 
limit. 
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We note that 



(15) 



with 



(16) 



Eq. (15) implies that we can solve the set of equations H^{fj,)rj{iJ,) = b for 
several masses, fi, simultaneously (for the same right hand b) using the multiple 
Krylov space solver described in Rcf. [15]. We will refer to this as the outer 
conjugate gradient inversion. It is easy to see that [i^^(/x),75] = 0, implying 
that one can work with the source b and solutions r]{ii) restricted to one chiral 
sector. 

The numerically expensive part of the Overlap-Dirac operator is the action 
of -f^o(O) on a vector since it involves the action of [y5e{Hyj) + e{Hyj)^5] on a 
vector. If the vector b is chiral (i.e. 756 = ±6) then, [y^€{Hyj) + e(i7^)75]6 = 
[75 ± l]e{H^)b. Therefore we only need to compute the action of e{Hyj) on a 
single vector. 

To study the possible onset of spontaneous chiral symmetry breaking in quenched 
QCD, we stochastically estimate, for a fixed gauge field background, 

l^'Z{i'{x)^{x))A = \^nb-\t,)] (17) 

and average over gauge fields. We also compute stochastically uj — x-k — Xao 
u; = |:(TV(75^)-^(//) + IV^-^(//)) = -(ViV')-2(-f (V'V')a) ■ (18) 

For a derivation of the above equation we refer the reader to Ref. [8]. Some 
simple manipulations yield the following relations 

{b\D-\n)\b)^:^^bHvin)-b) 
(6|(75^)-^(/^) + D-'{f^)\b) = ^^4^ (v'il^) - b') ivif^) - b) (19) 
where 

HoifJ^)vifJ^) = b with 756 = ±6 . (20) 
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Fig. 4. (a) CO without the contribution from global topology as a function of ji. (b) 
The pseudoscalar mass squared, arn?pg, and the vector mass, amy, for /3 = 5.85 and 
5.70 on an 8^ X 16 lattice as a function of /i. 

As discussed in Ref. [8], it is appropriate to remove the topological contribu- 
tions to the above quantities in order to study the onset of chiral symmetry 
breaking. For this we first compute the low lying spectrum of 75/^(0) using 
the Ritz functional method [22] . This gives us, in particular, information 
about the number of zero modes and their chirality. In gauge fields with zero 
modes we always find that all \Q\ zero modes have the same chirality. We 
have not found any accidental zero mode pairs with opposite chiralities. We 
then perform a stochastic estimate in the chiral sector that has no zero modes 
and double the result to get the total contribution to {ipip) and uj excluding 
topology. In this sector, the propagator is non-singular even for zero fermion 
mass. Given a Gaussian random source 6 with a definite chirality all we have 
to do is solve the equation H^(iJ,)rj{iJ,) — b for several values of 

In Fig. 4a we show u without the topology term added for various lattice 
sizes and (3 in SU(3) using m = 1.65. We see some indication of the onset of 
spontaneous chiral symmetry breaking (with strong finite volume dependence) 
at /3 = 5.85 where, as // decreases, there is a small region where a; ~ 1///, then 
ui turns over and goes like n^. This latter behavior is expected in finite volume 
and is obvious from the explicit /i^ dependence in Eq. (19). 

We show in Fig. 4b pseudoscalar and vector meson masses from a preliminary 
spectroscopy calculation for SU(3) /3 — 5.85 and 5.70 on an 8^ x 16 lattice. 

Masses are extracted using multiple correlation functions in an excited state 
fit. The fermion masses have been chosen to be above the region of decreasing 
LV from finite volume dependence in Fig. 4a, namely jj, > 10~^. As in the 
calculations for cu above, a multiple mass shift conjugate gradient solver was 
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Fig. 5. Plots of the probability distribution, Pmmiz), of the smallest (rescaled) eigen- 
value z for the orthogonal, unitary and symplectic ensembles in the lowest two 
topological sectors. Here the smallest eigenvalue A is rescaled by the volume V and 
the infinite volume chiral condensate T, to z = WS. The curve in each plot is a 
fit to the prediction from random matrix theory with the best value for the chiral 
condensate S. 

used for several values of fi in the solution of H^{ii)r]{fi) = b with chiral source 
b. We see some slight deviation of am^pg from linearity for decreasing /i. and 
arripg does not extrapolate to at = which we attribute to finite volume 
dependence. The vector mass my is fairly linear over the entire region. 



4 The Overlap-Dirac operator and random matrix theory 



The Goldstonc pions, associated with the spontaneous breaking of chiral sym- 
metry dominate the low-energy, finite-volume scaling behavior of the Dirac 
operator spectrum in the microscopic regime, defined by 1/Aqcd « L « 
with L the linear extent of the system. The properties in this regime 
are universal and can be characterized by chiral random matrix theory (RMT) 
within three ensembles, depending on some symmetry properties of the Dirac 
operator, and according to the sector of fixed topology, entering via the num- 
ber of exact zero modes (see [24] for a recent review). Since the Overlap-Dirac 
operator has the same chiral properties as the Dirac operator in the con- 
tinuum, and since it has exact zero modes in topologically non-trivial gauge 
fields, it is well suited to test the predictions of RMT. In Figure 5 the dis- 
tribution of the lowest (non-zero) eigenvalue is compared to the predictions 
of chiral RMT for examples in all three universality classes - SU(2) in the 



12 



fundamental representation for the orthogonal ensemble, SU(3) in the fun- 
damental representation for the unitary ensemble and SU(2) in the adjoint 
representation for the symplectic case - and in the sectors with zero or one 
exact zero modes [25]. Excellent agreement is seen. In addition, the condensate 
E obtained in the two different sectors of each ensemble from fits to the RMT 
predictions agreed within errors. This agreement further validates the chiral 
RMT predictions on the one hand and strengthens the case for the usefulness 
of the overlap regularization of massless fermions on the other hand. 



5 Conclusions 



The Overlap-Dirac operator provides a formulation of vector gauge theories on 
the lattice with an exact chiral symmetry in the massless limit and no fermion 
doubling problem. The use of the Overlap-Dirac operator is, however, CPU 
time intensive. We reviewed a few methods to implement the operator acting 
on a vector. Of these, we found the optimal rational approximation method, 
in conjunction with the exact treatment of a few low lying eigenvalues and 
eigenvectors of Hyj in e{Hyj), the most efficient. Further improvements in the 
numerical treatment of the Overlap-Dirac operator would be very helpful. 

The Overlap-Dirac operator has exact zero modes with definite chirality in 
the presence of topologically non-trivial gauge fields. Due to their good chiral 
properties overlap fermions are well suited for the study of spontaneous chiral 
symmetry breaking. It is possible to separate the contribution of the exact zero 
modes due to topology in a numerical computation. We have presented sample 
results in a quenched theory from the remaining non-topological modes. We 
presented first spectroscopy results with overlap fermions in quenched lattice 
QCD. Finally, we compared the distribution of the smallest eigenvalue of the 
Overlap-Dirac operator with the predictions from random matrix theory. 
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